Reconstructing the interaction rate in holographic models of dark energy 
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We reconstruct the interaction rate of the holographic dark energy model recently proposed by 
Zimdahl and Pavon in the redshift interval < z < 1.8 with observational data from supernovae 
type la, baryon acoustic oscillations, gas mass fraction in galaxy clusters, and the growth factor. It 
shows a reasonable behavior in the sense that it increases with expansion from a small or vanishing 
value in the long past but starts decreasing at recent times. The later feature suggests that the 
equation of state parameter of dark energy does not cross the phantom divide. 



I. INTRODUCTION 



Nowadays it is widely agreed that our Universe is currently accelerating its expansion -see [2| for recent reviews. 
However, this consensus does not extend to the agent behind this acceleration and in fact there are many 
competing candidates. The simplest one, a tiny cosmological constant, does nicely well at the pragmatic 
observational level- but it entails seemingly unsurmountable problems on the theoretical side. This is why many 
researchers are considering other possibilities, particularly some or other scalar of tachyon field with a strong 
negative pressure high enough to drive cosmic acceleration -these go under the collective name of "dark energy" 
fields. 

Among the most recent generic proposals there is a very suggestive one based on the holographic principle. 
Loosely speaking, the latter asserts that the entropy of a system is given by the number of degrees of freedom 
lying on the surface that bounds it, rather than in its volume The roots of this principle are to be found 
in the thermodynamics of black holes [I[ . Nevertheless, as noted by Cohen et al. 0] , a system may satisfy the 
holographic principle and, however, include states for which its Schwarzschild radius is larger than system size, 
L. This can be avoided by imposing the constraint that the energy of the system should not exceed that of black 
hole of the same size or, equivalently, p < 3 c 2 / (8ir G L 2 ), where c 2 is a (non-necessarily constant) parameter. 
In the cosmological context L is usually taken either as the event horizon radius or the Hubble radius. (For a 
quick summary of holographic dark energy see section 3 of Ref . [l| ) . 

The model of Ref. [1| is based in two main assumptions, (i) dark energy complies with the holographic principle 
with L identified as the radius of the Hubble horizon, ff -1 , hence p x — 3c 2 H 2 / (8nG) , and (ii) dark energy and 
dark matter do not evolve separately but they interact. Accordingly, the energy conservation equations are 

p m + 'iHp m = Q, p x + 3H(1 + w)p x = -Q , (1) 

where w is the equation of state parameter of dark energy, p x /p x , which is not constrained to be a constant. 
Subscripts m and x are for dark matter and dark energy, respectively. 

It should be noted that for spatially flat universes in the absence of interaction, Q = 0, there would be no 
acceleration [H, @]. Besides, Q must be a positive-definite quantity for the coincidenceproblem Q to be solved 
(or at least alleviated) [|J, and the second law of thermodynamics to be fulfilled 9]. Further, it has been 
forcefully argued that the Layzer-Irvine equation (lp| when applied to galaxy clusters reveals the existence of 
the interaction [ll|. To the best of our knowledge, the interaction hypothesis was first introduced, well ahead 
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of the discovery of late acceleration, by Wetterich [12j to reduce the theoretical huge value of the cosmological 
constant, and was first used in connection to holography by Horvat [HI]. As we write, the body of literature 
on the subject is steadily growing -see [l[ and references therein. Most cosmological models implicitly assume 
that matter and dark energy couple gravitationally only. However, unless there exists an underlying symmetry 
that would set Q to zero (such a symmetry is still to be discovered) there is no a priori reason to discard the 
interaction. Ultimately, observation will tell us whether the interaction exists. 

Following p}, we will write the interaction as Q = p^T, where T is an unknown, semi-positive definite, function 
that measures the rate at which energy is transferred from dark energy to dark matter. Clearly, as long as 
the nature of both dark ingredients of the cosmic substratum remain unknown, T cannot be derived from first 
principles; however, one can resort to observational data (in our case, supernovae type la (SN la), baryon 
acoustic oscillations (BAO), gas mass fraction in galaxy clusters and the growth factor) to roughly reconstruct 
it. The next section focuses on reconstructing the dimensionless quantity T/3H. 



II. RECONSTRUCTION 



The evolution equation 



(1 



3Hw- 



l + r 
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for the ratio r = p m /px between the energy densities follows from Eqs. |T]) and the above expressions for p x 
and Q. With the help of Friedmann equation Q m + Q x + Qk = U in terms of the usual density parameters 
Qi = 8TrGpi/(3H 2 ) (i — m,x), and f2fe = —k/(a 2 H 2 ), where k stands for the spatial curvature index of the 
Friedmann- Robertson- Walker metric, we can write 



Hi, 



(3) 



with q = —a/(aH 2 ), the deceleration parameter. Here we have assumed p x cx H 2 for holographic dark energy 
assuming the horizon as the Hubble horizon. 

Likewise, starting from the first of Eqs. |T]) and using Friedmann equation, we get for the equation of state 
parameter the expression 



w(z) = (1 + r) 
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where z denotes the redshift factor and a prime indicates derivative with respect to this quantity. 
We fit the Chevallier-Polarsky-Linder parametrization [14j , namely, 



w(z) — Wq + W\ 



1 + z' 



(•5) 



where wo is the present value of w(z), and wi a further constant, to current data from different observational 
probes and subsequently use the fitting values for wq and w\ to reconstruct the dimensionless ratio T/3H. 



As for the data, we resort to the various SN la observations in recent times. In particular we use 60 Essence 
supernovae 57 SNLS (Supernova Legacy Survey) and 45 nearby supernovae. We have also included the 
new data release of 30 SNe la detected by the Hubble Space Telescope and classified as the Gold sample by 
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Riess et al. 
192. 



[151 ]. The combined data set can be found in Ref. 16]. The total number of data points involved is 



Next we add the measurement of the CMB (Cosmic Microwave Background) acoustic scale at zbao = 0.35 as 
observed by the SDSS (Sloan Digital Sky Survey) for the large scale structure. This is the Baryon Acoustic 
Oscillation (BAO) peak) 



We also consider the gas mass fraction of galaxy cluster , f gas = M gas /M t ot, inferred from the X-ray observations 



[18| . This depends on the angular diameter distance (La to the cluster as f gas 
involved is 26. 



,3/2 



The number of data point 



Likewise, the 2dF galaxy redshift survey has measured the two point correlation function at an effective redshift 
of z s — 0.15. This correlation function is affected by systematic differences between redshift space and real space 
measurements due to the peculiar velocities of galaxies. Such distortions are expressed through the redshift 
distortion parameter j3. Correlation function can be used to measure it as j3 = 0.49 ± 0.09 at the effective 
redshift of z = 0.15 of the 2dF survey. This result can be combined with linear bias parameter b = 1.04 ± 0.11 
obtained from the skewness induced in the bispectrum of the 2dFGRS by linear biasing to find the growth factor 
g at z = 0.15, namely g = 0.51 ± 0.11 [Hj]. 



A. The spatially flat case 



The simplest case is when fifc = 0. Then, from (J3J) , r = ro , where the zero subscript means present value, and 
equations ([2]) and ([4]) reduce to 
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and 



w(z) = (1 + r Q ) 
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respectively. Using these two expressions we determine wq and w\ from the data and, with them, we reconstruct 
T/3H -see figures [T] and H 



The best fit values, with la error bars for the parameters when all the data (SN la + BAO + 
factor) are included, come to be: wo = —1.13 ± 0.24, w\ = 0.66 ± 1.35 (for fl m Q = 0.25 & flk 
w = -0.80 ± 0.28, wi = -1.75 ± 1.79 (for O m0 = 0.3 & fl k = 0, Fig. 



x-rays + growth 
= 0, Fig. [I]); and 



Here, one cautionary remark seems in order. The fact that r was never large might lead the reader to think 
that the model of Ref. [l| seriously conflicts with the standard scenario of cosmic structure formation. One 
may believe that at early times the amount of dark matter would have been too short to produce gravitational 
potential wells deep enough to lead to the condensation of galaxies. However, this is not so; a matter dominated 
phase is naturally recovered since at high and moderate redshifts the interaction is even smaller than at present 
whence the equation of state of the dark energy becomes close to that of non-relativistic matter -see [l[ for 
details. 



A related point is to realize that dark energy clusters similarly to dark matter when the equation of 
state of the former stays close to that of the latter. In this connection, it is worthwhile to recall the 
perturbation dynamics of this model. This was studied in 1] making use of the perturbed metric 
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FIG. 1: The dimensionless ratio V/(3H) vs redshift. In the four panels we have fixed Q,„o = 0.25 and Q,k = 0. The solid 
line is for the mean value and the shaded area indicates the la region. The region above the horizontal dashed line can 
be visited only when the dark energy becomes of phantom type, i.e., w < — 1. 

ds 2 = — (1 + 2tp) dt 2 + a 2 (1 — 2ip) S a f3dx a dx^ , with -0 the scalar metric perturbation, and the Bardeen 
gauge-invariant variable [20] 



C = -i> + 
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which represents curvature perturbations on hypersurfaces of constant energy density. Here, an upper-hat means 
perturbation of the corresponding quantity; likewise, p = p m + p x and p = p x . 

Corresponding quantities for the components are 



-tj) - H— 

PA 



(A = m, x). 



(9) 



On large perturbation scales we have that 



(10) 



with P = p/(p + p) and D = p/(p + p) and parallel expressions for Cm and ( x . Therefore, insofar as both 
equations of state do not differ significantly the evolution of these two perturbations will be alike. 



In the particular, simplest, case of T = constant last equation reduces to 
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FIG. 2: Same as Fig. 1 except that here we have fixed Q m o = 0.30 and Q,k = 0. 



s 6(l + £)r 2 ' 

This one can be readily integrated to (see section 6 of Ref. [l| for details) 



(11) 



c = c, 



Tr 



1 



3 r 3H,r 
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Again, so long as the the equation of state parameter of dark energy w remains close to that of dark matter, both 
components will cluster in a similar fashion. Further, a non-vanishing interaction introduces a non-adiabatic 
feature that grows as a 3 / 2 which will have an impact on the integrated Sachs- Wolfe effect. Possibly this feature 
might be used in the future to discriminate the model under consideration from the ACDM model -recall that 
in the latter £ remains constant and does not produce the said effect. 



B. Non-spatially flat cases 



When Qk =/= the ratio r between energy densities is no longer a constant. This is an extra unknown function 
in our fitting procedure. But one should not expect a large variation in r in the redshift range that has 
been considered in this paper. In our subsequent computation, we take the Taylor series expansion for r 
around its present day value and take up to the first order term in the expansion. We therefore parameterize it as 

2: 

r = r + ri(l-a)=r + ri— — , (13) 
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FIG. 3: Same as Fig. 1 except that we have fixed Q m o = 0.30 and fifco = 0.002. 



with tq is the present day value for r. r\ is a constant which can be related to the present ratio of densities 
between flk and Q x by 



n 

2 



IT 
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This can be used to fix the unknown constant r\ for a given flko and fl x a- Also 



1 ( ,1 + z 

r — wr 



3H 1 + r 
where w is given by Eq. ([?]). 

Using these expressions, the ratio T/3H is reconstructed from the data in Fig. [3J 



(15) 



The best fit values, with lcr error bars for the parameters when all the data (SN la + BAO + x-rays + growth 
factor) are included, come to be: wo = —0.806 ± 0.29, w\ = —1.74 ± 3.33. It is seen that the curvature, being 
small as WMAP 3yr [2l| tells us, has little consequence on the evolution of the interaction rate. 
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III. CONCLUDING REMARKS 



We reconstructed the interaction term Q of Ref. [l| in the redshift interval (0 < z < 1.8) of observational data 
(supernovae type la, baryon acoustic oscillations, gas mass fraction, and growth factor). The interaction rate 
T (and hence Q) is always positive, its general trend is to decrease as z increases but it shows no indication of 
becoming negative at larger redshifts. This corroborates that as previously suggested 0, [Hj] the energy transfer 
proceeds from dark energy to dark matter rather than otherwise. While phantom behavior cannot be excluded 
at recent and present times it only occurs in a manifest way either for large fi x o -see Fig. [TJ- or when just the 
supernovae data are considered (top-left panel of Figs. [US])- When fl x0 is a bit lower (say, 0.7) and BAO and 
other data are included, the mean value of dimensionless interaction rate, T/3H, no longer crosses the phantom 
divide (i.e., the horizontal dashed line). It simply reaches a maximum near z = and decreases with expansion. 
This is a reassuring result as holography is not compatible with phantom energy [23]. On the other hand, it 
should be noted that O x o values as high as 0.75 do not seem favored from a combination of results from WMAP 
lyr and weak lensing which yields Q x q — 0.70 ± 0.3 (23j. 

Adding a small curvature term -say, ttko = 0.002-, has only a tiny impact (compare Figs. [2] and [3j. This also 
holds true when the curvature bears the opposite sign; this is why we have not included a corresponding figure. 

In any case, it should be noted that the concordance ACDM model (wo = —1, Wi = 0) shows compatibility 
within ler confidence level with the set of data considered in this work. 
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